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We calculate the Hall conductivity a xy and resistivity p xy of a granular system at large tunneling 
conductance jt > 1. We show that in the absence of Coulomb interaction the Hall resistivity 
depends neither on the tunneling conductance nor on the intragrain disorder and is given by the 
classical formula p xy = H/(n*ec), where n* differs from the carrier density n inside the grains by a 
numerical coefficient determined by the shape of the grains. The Coulomb interaction gives rise to 
logarithmic in temperature T correction to p xy in the range T < T < Txan(grE c , Eth): where V is 
the tunneling escape rate, E c is the charging energy and £rh is the Thouless energy of the grain. 

PACS numbers: 73.63.-b, 73.23.Hk, 61.46.Df 



Hall resistivity (HR) of metals and semiconductors 
gives very important information about their properties. 
According to the classical Drude-Boltzmann theory HR 
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H/(nec) 
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does not depend on the mean free path and allows one 
to experimentally determine the carrier concentration n. 

Recently, much attention from both experimental and 
theoretical sides has been paid to granular systems (see 
a review []J and references therein). Although vari- 
ous physical quantities have been calculated in different 
regimes, the Hall transport in the granular matter has 
not been addressed theoretically yet. In this work we 
calculate the Hall conductivity (HC) of a granular system 
and Coulomb interaction corrections to it in the metal- 
lic regime, when the intergrain tunneling conductance 
Gt = (2e 2 /h)gT is large, gr ^> 1 (further we set K = 1). 

Technically, calculating HC a xy appears to be more 
complicated than calculating the longitudinal conductiv- 
ity (LC) a xx . The granularity of the system is ensured 
by the condition that the conductance Go = 2e 2 <?o of the 
grain is much larger than the tunneling conductance Gt, 
go/gr 3> 1. In this limit the main contribution to o~ xx 
comes from the tunnel barriers between the grains rather 
than from scattering on impurities inside the grains. In 
the absence of Coulomb interaction LC equals 



^ = G T a 2 -«, 



(2) 



where a is the size of the grains and d is the dimensional- 
ity of the system. Therefore when studying longitudinal 
transport one can neglect electron dynamics inside the 
grains, which simplifies calculations significantly. On the 
contrary, for Hall transport one is forced to take the in- 
tragrain electron dynamics into account, since the Hall 
current originates from the transversal drift in crossed 
magnetic and electric fields inside the grains. 

As we find in this work the intragrain electron dynam- 
ics can be included within the diagrammatic approach by 
considering higher diffusion modes inside the grain. This 



procedures accounts for the finitcncss of the ratio gr/ga, 
and allows one, in principle, to study both LC and HC 
of the granular system for arbitrary ratio gx / 'go- The ob- 
tained results reproduce the solution of the classical elec- 
trodynamics problem for a granular medium (e.g. the 

formula ai°J — a 2 ~ d GTGo/ {Gt + Go) can be obtained 
as a series in <?t/.9o)- Quantum effects (e.g. Coulomb 
interaction and weak localization) may be incorporated 
into this scheme afterwards. 

We perform calculations for magnetic fields H such 
that wjt <C 1, where ujj = eH/mc is the Larmor fre- 
quency and r is the electron scattering time inside the 
grain. Since the effective mean free path I = VpT < a 
does not exceed the grain size a, and typically a « 
10 — lOOnm, the condition ujt <§; 1 is well fulfilled for 
all experimentally available fields H. 

Let us list the main results of this work. First we 
neglect Coulomb interaction completely and obtain for 
HC in the lowest nonvanishing order in gr/go 



a£°) = G 2 T R H a 
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where Rjj is the classical Hall resistance of a single grain 
(see Fig. l(right)). This result obtained by diagrammatic 
methods is completely classical provided the tunneling 
contact is viewed as a surface resistor with conductance 
G T . The HR of the system 



Pxy = RhCL 
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= H/(n* d ec) 



(4) 



is given by the Hall resistance of a single grain Rjj , which 
depends on the geometry of the grain but not on the intra- 
grain disorder. Eq. defines the effective carrier den- 
sity n* d of the granular medium. For a three-dimensional 
(d = 3, many granular monolayers) array = An differs 
from the electron density n in the grain by a numerical 
factor A, < A < 1, determined by the grain geometry 
[l4| . For grains of a simple geometry this factor is given 
by the ratio of the largest cross section area S to the cross 
section area of the lattice cell a 2 : A = S/a 2 . So, A = 1 
for cubic grains (S — a 2 ), A = ir/4 for spherical grains 
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(S = 7ra 2 /4). For a two-dimensional (d = 2, single gran- 
ular monolayer) array the 3D result must be multiplied 



by the thickness of the layer a: 



= aAr 



Next we calculated the first-order correction to HC 
a x ° y (Eq. ©) due to Coulomb interaction at tempera- 
tures T > r not much smaller than the tunneling escape 
rate T = gxS (6 is the mean level spacing of the grain). 
We find significant T-dependent corrections in the range 
r < T < QtE c (E c = e 2 /a is the charging energy of the 
grain), whereas for T > grE c the relative corrections are 
of the order of 1/gr or smaller . 

There exist two different contributions to HC 
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One of them, Sag , can be attributed to the renormaliza- 
tion of individual tunneling conductances Gt (tunneling 
anomaly 0,0,0]) in a granular medium and has the form: 
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FIG. 1: (left) Diagrams for the Hall conductivity al^J of the 
granular system (Eqs. 1141 1. The diffuson Dy connecting 
contact Si to so is shown. External tunneling vertices (wavy 
lines) must be attached in 4 possible ways. Three more dif- 
fusons D^, Dy , connecting contacts S2,S3,S4 to so, re- 
spectively, must be also taken into account, (right) Classical 
picture of Hall conductivity of a granular system. The current 
I y = GrVy running through the grain in y-direction causes 
the Hall voltage drop Vh = RhI v between its opposite banks 
in x-direction, which is also applied to the contacts (the total 
voltage drop per lattice period in i-direction is 0) giving the 
Hall current I x = G T V H = G 2 T R H V y (Eq. @). 



(2) 

The other one, Sa xy , is due to elastic cotunneling 
(EC) it involves diffusion of an electron through the 
grain (that is why it is suppressed at temperatures 
greater than the Thouless energy of the grain £rh): 
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where Cd is a numerical factor Eq. i|16|) . Since the cor- 
rection Sag merely renormalizes the conductance Gt in 
Eq. ©, it does not affect the HR p xy = a xy /a xx . Indeed 
the correction Eq. JBJ is cancelled by the corresponding 
correction to a xx , describing renormalization of 
Gt in Eq. (J2J) • Therefore the total correction Sp xy to HR 

P X y — Pxy 

+ Sp xy is due to the EC effect (Eq. 0) only: 
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(2) 



o (0) 

Pxy 



7 (0) 
>xy 



Cd 
4irg T 



In 



mm(g T E c , E Th ) 



T 



(8) 



We conclude that the Hall resistivity p xy of a granu- 
lar metal at temperatures T > miH-(g r pE cl -£/xh) is given 
by Eq. (0} and is independent of the intragrain and tun- 
nel contact disorder. Measuring p xy at such T and using 
Eq. 0} one can extract an important characteristic of the 
granular system: its effective carrier density n* d . At tem- 
peratures r < T < miii(gTE c , i/rh) Coulomb interaction 
leads to In T-dependent corrections to p xy . Comparison 
of Eqs. (QJ-JSJ) with experimental data may serve as a 
good check of the theory developed here. The logarith- 
mic dependence p xx = a + b In T of granular metals has 
been observed experimentally ||, and p xy can also be 
measured [t| . Our theory may also be applied to indium 
tin oxide(ITO) materials (see e.g. ^(J). Another related 



effect is the anomalous Hall effect in ferromagnetic gran- 
ular materials [ll| . 

Now we briefly outline the model and method used 
to derive the announced results, details of our calcu- 
lations will be presented elsewhere |l2j. We consider a 
quadratic(d = 2) or cubic(d = 3) lattice of identical (in 
form and size) grains coupled to each other by tunnel con- 
tacts (Fig. 1). To simplify the calculations we assume the 
intragrain electron dynamics diffusive, I <C a. However, 
our results are also valid for ballistic (I ~ a) intragrain 
disorder. In the metallic regime (gr ^S> 1) the Coulomb 
interaction effects can be considered as a perturbation. 

We write the Hamiltonian describing the system as 

H = H o + H t + H c . 



Hn = 



]T J dr^rOfc (pi-e/cA(r,)) + tf(r,)]^(r,) 



is the electron Hamiltonian of isolated grains, £(p) = 
p 2 /2m — cf, A(ri) is the vector potential describing uni- 
form magnetic field H = He z , U(ri) is the random dis- 
order potential of the grains, i = (ij, . . . , id) is an integer 
vector numerating the grains. The disorder average is 
performed using a Gaussian distribution with the vari- 
ance (J/(ri)[/(rQ) = 2 1 8{ti — r-), where v is the den- 
sity of states in the grain at the Fermi level per one spin 
direction. Further, the tunneling Hamiltonian is given by 



(iJ> 



E 



dsidsji(si,Sj)V> T (si)V>(Sj) 
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the summation is taken over the neighboring grains con- 
nected by a tunnel contact, and the integration is done 
over the surfaces of the contact. The inevitable irregu- 
larities of the tunnel contact are modelled by the relation 
(t(si,Sj)t(sj,Si)) = ig<5(si — sj), where 5(sj — sj) is a 6- 
function on the contact surface, has the meaning of 
tunneling probability per unit area of the contact. The 
Coulomb interaction can be taken in the form (e.g. 

fr e2 V- m»J 



where m — J dri t/;' (ri)ip(r{) — n is the excess number of 
electrons in the i-th grain, Cy is the capacitance matrix. 
The conductivity is calculated using the Kubo formula: 



(uj) = 2eV- d i^(n ab (w,i-j)-n ab (0,i-j)) (9) 



Cab 



where u> = 2irTk > is a bosonic Matsubara frequency, 
a and b are the lattice unit vectors, 

n ab (o;,i-j) = f 1 dre^ r (r T / i!a (T)/ jib (0)) (10) 
Jo 

is the current-current correlator, /^(t) = Ai +a ,i(r) — 
Xi,i +a (r), the thermodynamic average (. . .) is taken with 
H, A(t) = e Hr Ae~ Hr is the Heisenberg operator. 

We perform calculations of <r a b using diagrammatic 
technique. Let us first neglect Coulomb interaction H c . 
Technically, one expands Eq. (|10fl both in disorder po- 
tential U(ri) and the tunnelling Hamiltonian H t and av- 
erages over the intragrain and contact disorder. 

Summing the diagrams corresponding to intragrain 
motion we obtain an important object: the diffusion 
propagator (diffuson) of a single isolated grain given by 
a disorder-averaged product of two Green functions: 

£>(w,r,r') = — (G(e+w,T,r')G(e,r',r)), (s+uj)e < 0. 

I'KV 

In the presence of magnetic field (luht <C 1) this propa- 
gator satisfies the equation 

(M-DoV?)D(w,r,0 = J(r-r') (11) 

and the boundary condition at the grain surface 

(n, V r D)\ reS = WffT(t, V r D)\ reS - (12) 

Here Dq = vpl/3 is the diffusion coefficient in the grain, 
n is the normal unit vector pointing outside the grain, 
t = [H, n]/|[H, n]| is the tangent unit vector pointing 
in the direction of the edge drift. Eq. (|12|l is due to 
the fact that the current component normal to the grain 
surface vanishes, the right hand side describes the edge 



drift caused by magnetic part of the Lorentz force. The 
solution to Eqs. i|ll|) . l|12[l can be presented in the form 



^(rR(r') 



(13) 



where </>„ are the eigenf unctions of the problem 

-D a V 2 r 4> n = j n 0n, (n, V r </>„)|s = w ff r(t, V r (/> n )|s- 



There always exists a uniform solution 4>o(r) = 1/W (V 
is the grain volume) with zero eigenvalue 70 = giving 
the 0-th harmonic l/|w|V in Eq. 

Accounting for tunneling processes makes the trans- 
port between different grains possible. The diagrams for 
Eq. JTJJ in the 2iV-th order in H t are given by diffusons 
connecting the contacts (i + a, i) and (j + b, j) that run 
through the other N — 1 contacts in all possible ways. 
An important observation is that the diffusons inside the 
grain always enter the expression for cr a b as a difference 
D(oj, Si, S2) — D(lj, S3, S4) of the diffusons connecting dif- 
ferent contacts. Therefore the 0-th harmonic drops out 
and the contribution to er a b comes from nonzero harmon- 
ics, each 2 powers of Ht bringing a factor T / E^h = gx /go- 

For the longitudinal conductivity (a = b = e^) the 
above procedure would lead to a small contribution in 
the limit gr/go 1, since Eq. J2Jl can be obtained from 
a single contact (j = i) without expanding Eq. (|10f> in 
H t (N — 0). For the Hall conductivity (a = e^b = 
e^) considering the higher diffusion modes is inevitable, 
because even in the lowest nonvanishing order in gr/go 
we have to connect different contacts by diffusons of a 
single grain as shown in Fig. l(left) (j = i, i + e x , i — e y , 
i + e x — e y ). Doing so, we obtain 



4°» 
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D / -Dk 



(14) 



where gx — 2Tr(vtQ) 2 Y l is the conductance of a tun- 
nel contact, S is the area of the contact, D a — 
^2 J ds ds a D(s ,s a ), a = 1,2,3,4 for a =/\\,^,\ 
respectively (Fig. l(left)). Here 

D(r,r') = ^^(r)^(r')/7n 

n>0 

is the diffuson without the 0-th harmonic at u> = 
satisfying Eqs. lfTT|) . l(T^|l with u — 0. Retaining 
only the 0-th harmonic in Eq. (|13|l would give just 
in Eq. 114JI and we are forced to take all higher har- 
monics into account. Eq. (|14|l is nonzero for H 7^ 
since the edge trajectories for Dy = Dy are shorter (if 
e > is assumed) than for = D\, and therefore 
Dy — D\ = Dy - D\ > 0. 

In fact, the result Eq. I|14|l is purely classical, provided 
one treats the tunnel contact as a surface resistor with 
the conductance Gt- Indeed, the classical HC of the 
granular medium can be easily presented in the form of 
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Eq. © (Fig. 2 (right)). The Hall resistance R H of the 
grain is defined via the difference (Hall voltage) of elec- 
tric potential <p(r) between the opposite banks of the 
grain, Vh = ¥?(s r ) — ¥>( s z) = Rtily, when the current 
I y = I passes through the grain. The current density 
j(r) = — <7 V<p(r) (cto is the conductivity tensor) satis- 
fies the continuity equation divj = q(r) and the bound- 
ary condition (n • j)|s = 0. The charge source func- 
tion q(r) is nonzero on the contacts surface only and 
/ ds d q(s d ) = — J ds u q(s u ) = I. Therefore ip(r) satisfies 

-V 3 t ip = q(r)/af x , (n, V<p)\ 8 = Wjrr(t, V<p)\ s . (15) 

Comparing Eq. (|T5|) with Eqs. l(TT )l .l(T^ )l we see that 
D(r, r') is a Green function for the problem Eq. Q15[l. 
Thus the solution to Eq. (|15|l is 



(p(r) = 



1 I 



ds d D(r,s d )- / ds u -D(r,s 



and Eq. (JTIJl leads to Eq. © (Einstein relation crf^. = 
2e 2 vD§ was used). This establishes the connection be- 
tween our diagrammatic approach of considering higher 
diffusion modes and the solution of the classical electro- 
dynamics problem for the granular system. 

Luckily, for simple geometries (cubic, spherical) of the 
grain the Hall resistance Rh can be obtained from sym- 
metry arguments without solving the problem Eq. (|15fl . 
In such cases the Hall voltage equals Vh = pfLal/S, 
where S is the area of the largest cross section of the 
grain, pf = H/ (nee) is the specific HR of the grain ma- 
terial expressed in terms of the carrier density n inside 
the grain. Therefore, Rh = pf y a/S and the HR of the 
granular medium can be expressed in the form of Eq. 
where n* d = a d ~ 3 An, A = S/a 2 < 1. The quantity n* d de- 
fines the effective carrier density of the system. 

Diagrammatic approach allows us to incorporate quan- 
tum effects of Coulomb interaction on Hall conductivity 
into the developed scheme. We omit details here leaving 
them for a more comprehensive version [l2^| . We consider 
the range of not very low temperatures T > T where 
we can neglect the large scale contributions analogous 
to those for homogeneously disordered metals [§]• We 
find that in this regime in the first order in the screened 
Coulomb interaction two contributions exist (i = 1,2): 



2e 2 a 



n>0 



where f n>a = ^ J ds ds a (p n (so)<f>n (s a ), a = 1,2,3,4 for 
a =y,\, ,/ ,\, respectively (Fig. 1 (left)) , and 



2 E 

n 



v Q (n) + v 2 {n)-2v 1 (n) 
{uj + n + j n )n 2 



Here 27rT£ n '.F(fi) = Eo<n<* SIF(Q) + E„<n^l 
V = V(i, i), Vi = V{i + e x , i), V 2 = V(i + e x + e y , i) 
are the components of the screened Coulomb interaction 

v(n,i,j) = j ^gj e Wi-j)v(n )q ) with 

V((l, q) - E c (q)/[l + 2E c (q)Tj(S\n\)}, 

T q = 2rE /3 (l - cos^o), E c (q) = e-^-fieVd-j 
(/? = x, y for d = 2 and f3 = x,y, z for d = 3, f2 = 

2TrTm, m G Z). The contribution Sa^J renormalizes 
individual tunneling conductances Gt in Eq. © , whereas 

(2) 

the contribution 5ox y is due to elastic cotunneling (note 
the diffuson + O + 7 n ) of the grain in the expression 
for Xn^)- The large logarithms in A« (u>) and \n\ui) 
come from frequencies < grE c , for which the Coulomb 
potential is completely screened: V(f2,q) = <5|£l|/2r q . 
Performing analytical continuation and extracting Oxy 
with the help of Eq. fH}l we arrive at Eqs. ©,0 with 



a d q (1 — cosqj;a)(l — cos q y a) 



[2n) a 



Y,p0- - cosqpa) 



(16) 



In conclusion, we presented a theory of the Hall con- 
ductivity of granular metals. We have shown that at 
high enough temperatures the Hall resistivity is given by 
a classical expression, from which one can extract the 
effective carrier density of the system. At lower tem- 
peratures charging effects give a logaritmic temperature 
dependence of Hall resistivity. We hope that our predic- 
tions may be rather easily checked experimentally. 

Authors thank I. S. Beloborodov and A. F. Volkov for 
illuminating discussions and acknowledge financial sup- 
port of Degussa AG (Germany), SFB Transregio 12, the 
state of North-Rhine Westfalia, and European Union. 



■/n,\+/n,^-/n,\)A«(a;) 



[1] 

[2] 

[3] 
[4] 
[5] 

[6] 

[7] 

[8] 

[9] 



I. S. Beloborodov, K. B. Efetov, A. V. Lopatin, and V. 
M. Vinokur, cond-mat/0603522 (2006). 
B.L. Altshuler and A.G. Aronov, in Electron- Electron In- 
teractions in Disordered Systems, edited by A.L. Efros 
and M. Pollak (Elsevier, Amsterdam, 1985). 
Yu. V. Nazarov, Sov. Phys. JETP Lett. 49, 126 (1989). 
L. S. Levitov, A. V. Shytov, JETP Lett. 66, 214 (1997). 
D. V. Averin and Yu. V. Nazarov, Phys. Rev. Lett. 65, 
2446 (1990). 

K. B. Efetov and A. Tschersich, Europhys. Lett. 59, 114, 

(2002); Phys. Rev. B 67, 174205, (2003). 

I. S. Beloborodov, K. B. Efetov, A. V. Lopatin, and V. 

M. Vinokur, Phys. Rev. Lett 91, 246801, (2003). 

A. Gerber, A. Milner, G. Deutscher, M. Karpovsky, and 

A. Gladkikh, Phys. Rev. Lett. 78, 4277, (1997); Simon, 

R. W., B. J. Dalrymple, D. V. Vechten, W. W. Fuller, 

and S. A. Wolf, Phys. Rev. B 36, 1962, (1987) 

X. X. Zhang, Chuncheng Wan, H. Liu, Z. Q. Li, Ping 

Sheng, J. J. Lin, Phys. Rev. Lett. 86, 5562 (2001); Peng 

Xiong, Gang Xiao, J. Q. Wang, John Q. Xiao, J. Samuel 

Jiang, C. L. Chien Phys. Rev. Lett. 69, 3220 (1992). 



5 



[10] J. Ederth et al, Phys. Rev. B 68, 155410 (2003); J. Ederth 

et al, Thin Solid Films. 445, 199 (2003). 
[Ill P. Mitra, A. F. H ebard, K. A. Muttalib, P. Woelfle, 

|cond-mat /0606215 
[12] M. Yu. Kharitonov, K. B. Efetov, in preparation. 
[13] I. L. Aleiner, P. W. Brouwer, and L. I. Glazman, Phys. 

Rep. 358, 309 (2002). 
[14] Note that although the granular array may be two- (d — 

2) or three-dimensional (d = 3) , the grains themselves are 



three-dimensional, and n is a three-dimensional density. 

[15] We do not consider weak localization effects in this paper 
assuming they are suppressed by either magnetic field or 
inelastic processes. 

[16] To get a correct aj-dependence of Eq. I14H (its indepen- 
dence of w) one should renormalize the diffusons by the 
Coulomb interaction inside the grain. This results in the 
replacement l/(|w| +7n) — > l/-y n 0- 



